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other  two  phases  are  written  in  MPROLOG  and  run  on  a  Motorola  68020-based 
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The  lists  of  bump  descriptions  showed  the  key  differences  between  two  different  ships  of 
the  same  ship-class,  indicating  future  programs  could  identify  ships  using  the  output  from 
the  programs  of  this  thesis. 


j  Mdfo 
ACKNOWLEDGEMENTS 


I  wish  to  express  my  deepest  appreciation  to  my  wife,  Sally,  for  her  encouragement 
and  support  during  the  research  and  writing  of  this  thesis.  Thank-you  also  to  my  three 
daughters  for  their  understanding  that  "daddy  had  to  go  to  school  a  lot." 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 1 

A.  OBJECT  RECOGNITION  AND  IDENTIHCATION 1 

B.  MILITARY  REQUIREMENTS  FOR  SHIP  IDENTIFICATION 1 

C.  DESCRIPTION  OF  THESIS 2 

II.  RESEARCH  IN  COMPUTERIZED  OBJECT  IDENTIHCATION 3 

A.  BACKGROUND 3 

B.  RESEARCH  USING  CONVENTIONAL  METHODS 3 

1 .  Shape  Recognition  Using  Binary  Images 3 

2.  Panem  Classification  from  Elemental  Shapes 4 

3 .  Identification  of  Spherical  Objects 4 

4 .  Parallel  Processors  and  Image  Analysis 5 

C.  RESEARCH  USING  AI  TECHNIQUES 5 

1 .  Knowledge-Based  Airport  Identification 5 

2.  Rule-Based  Methods  for  Image  Analysis 6 

3.  Recognizing  Man-Made  Objects  in  Aerial  Images 7 

III.  METHODS  OF  SHIP  IDENTIHCATION 8 

IV.  A  PICTURE-DESCRIPTOR  EXTRACTION  PROGRAM 10 

A.  OVERVIEW 10 

B.  FINDING  THE  SILHOUETTE  BOUNDARY 11 

C.  LOCATING  THE  BUMPS 13 

D.  BUMP  FEATURE  IDENTIFICATION 15 

1.  Removing  the  Shape  from  the  Silhouette 15 

2 .  Shape  Analysis  and  Identification 18 

E.  FINAL  OUTPUT  OF  THE  PROGRAM 18 

V.  RESULTS  OF  PICTURE-DESCRIPTOR  EXTRACTION  PROGRAM 19 

A.  FINDING  THE  SILHOUETTE  BOUNDARY 19 

B.  LOCATING  THE  BUMPS 20 

C.  BUMP  FEATURE  IDENTIFICATION 21 

D.  PROGRAM  LIMITATIONS  DUE  TO  HARDWARE  CONSTRAINTS 25 


VI.    CONCLUSION 27 

A.  SUCCESS  OF  FEATURE  EXTRACTION  AND  IDENTIFICATION 27 

B.  AREAS    FOR    FUTURE    WORK 27 

APPENDIX  A  -  C  PROGRAM  SOURCE  CODE 29 

APPENDIX  B  -  MPROLOG  PROGRAM  SOURCE  CODE 36 

APPENDIX  C  -  OUTPUT  FROM  GARCIA  (LAMPS) 55 

LIST  OF  REFERENCES 58 

INITIAL  DISTRIBUTION  LIST 59 


I.    INTRODUCTION 

A.  OBJECT  RECOGNITION  AND  IDENTIFICATION 

The  process  of  visually  identifying  an  object  can  be  broken  into  detection,  delineation, 
and  identification.  Detection  locates  the  object,  delineation  separates  the  region  the  object 
occupies  from  the  background,  and  identification  determines  what  the  object  is.  The 
identification  process  is  of  critical  importance  in  military  applications.  Human  analysts 
remain  the  only  means  of  performing  this  task  accurately  and  reliably. 

Computers  can  provide  assistance  in  each  of  these  three  areas,  although  they  are 
weakest  in  identification.  Their  processing  speed  allows  for  rapid  scanning  of  an  area  for 
detection  targets.  Once  detected,  there  are  many  ways  of  extracting  the  target's  shape. 

B.  MILITARY  REQUIREMENTS  FOR  SHIP  IDENTIFICATION 

Rapid  detection,  analysis,  and  identification  of  ships  is  important  in  military 
operations.  Identification  is  necessary  to  determine  whether  the  ship  is  friendly  or 
potentially  hostile.  The  number  of  action  options  available  to  a  commander  rapidly 
decreases  the  longer  it  takes  to  identify  the  target. 

Training  humans  in  ship  identification  is  time-consuming  and  expensive.  Once 
trained,  these  skills  suffer  if  not  practiced  frequently.  Additionally,  the  amount  of 
information  to  be  interpreted  can  be  overwhelming  and  the  time  constraints  minimal. 
Computer- assisted  identification  promises  to  alleviate  these  current  shortcomings. 

Our  approach  to  computer-assisted  ship  identification  is  to  break  a  detected  object 
down  into  its  elementary  pieces  and  then  identify  the  pieces.  A  digitized  ship  silhouette's 
boundary  is  delineated  by  a  threshold  comparison.  The  boundary  is  converted  into  a  list  of 
line-segment  coordinates  of  contour  turn  points.  Using  artificial-intelligence  techniques, 


these  points  are  sequentially  examined  to  find  and  extract  the  coherent  pieces  of  the  ship. 
The  size,  shape,  amount  of  irregularity,  and  location  of  each  extracted  piece  are  then 
examined.  Each  piece  is  given  an  identity  of  gun-turret/weapons-system,  superstructure, 
weapons-system,  mast/support,  antenna,  radar,  or  unknown.  The  list  of  these  pieces  and 
their  identity  can  be  reported  to  humans  or  could  be  supplied  to  future  systems  similar  to 
[Ref.  1]  which  could  identify  the  ship  from  its  known  pieces  and  a  database  search. 

C.    DESCRIPTION  OF  THESIS 

The  remaining  chapters  of  this  thesis  present  background  information  and  describe  the 
identification  programs  and  their  results.  Chapter  II  describes  some  of  the  current  and  past 
research  in  computer-assisted  imagery  and  object  identification.  Chapter  III  examines  the 
problems  inherent  in  object  identification  and  describes  areas  where  computer-assisted 
systems  would  improve  processing  speed  and  accuracy. 

Chapter  IV  is  a  discussion  of  the  programs  we  developed  for  feature  extraction  and 
identification,  the  algorithms  we  used,  and  the  output  the  final  program  generates.  Chapter 
V  will  examine  the  results  of  the  program  for  different  ship  images,  and  chapter  VI  presents 
the  conclusions  and  areas  for  follow-on  research. 


II.  RESEARCH  IN  COMPUTERIZED  OBJECT 
IDENTIFICATION 

A.  BACKGROUND 

Computer-automated  and  computer-assisted  object  recognition  and  identification  has 
applications  in  many  different  areas,  including  robotics,  imagery  interpretation,  and 
chemistry  [Ref.  2:  p.ll].  For  example,  imagery  analysts  need  help  in  processing  the 
large  quantity  of  digital  imagery  data  they  receive.  Much  of  their  work  involves  detecting 
changes  in  images  taken  at  different  times,  a  time-consuming  task  well-suited  to 
computerized  scanning  [Ref.  3:p.  905].  Programs  which  enable  computers  to  identify 
objects  have  involved  conventional  procedural  image  analysis,  artificial  intelligence  (AT),  or 
a  combination  of  both. 

B  .   RESEARCH  USING  CONVENTIONAL  METHODS 

1 .  Shape  Recognition  Using  Binary  Images 

Grogan  [Ref.  4]  performed  a  comparative  study  of  shape  recognition  and 
description  using  binary  images.  His  test  objects  were  six  different  aircraft  arranged  in 
different  views.  His  analysis  concentrated  on  the  use  of  global  shape  methods  to  identify 
and  analyze  the  image  boundary.  The  five  methods  compared  were  1)  Fourier  descriptors 
of  the  boundary,  2)  Walsh  points  of  the  boundary,  3)  the  cumulative  angular  deviant 
Fourier  descriptors,  4)  moments  of  the  silhouette,  and  5)  moments  of  the  boundary. 

The  research  was  conducted  in  1982  and  1983,  and  the  speed  of  the  programs  was 
good.  One  of  their  limitations  was  the  resolution  capabilities  of  the  equipment;  images 
were  limited  to  256  by  256  pixels,  so  silhouettes  were  simplified  drawings  of  the  aircraft  to 
be  identified. 


2.  Pattern  Classification  from  Elemental  Shapes 

Research  which  involved  dissecting  silhouetted  pictorial  patterns  into  specific 
elemental  shapes  was  conducted  by  Todd  [Ref.  5].  His  efforts  concentrated  on  pattern 
classification.  A  figure  was  dissected  into  subfigures  and  each  subfigure  was  analyzed 
according  to  multiplicity,  orientation,  position,  size,  shape  and  position.  The  analysis  of 
the  subfigures  generated  a  Figure  Classification  Number  (FCN)  for  the  whole  pattern, 
which  could  be  compared  against  other  known  FCN's  in  a  database  search.  A  match 
would  allow  the  figure  to  be  identified. 

The  subfigures  of  an  object  generated  by  Todd  were  not  each  classified  as  a  named 
feature  of  the  object,  but  the  analysis  of  these  pieces  led  to  an  object  identification.  The 
program  successfully  recognized  22  different  aircraft  silhouettes  as  airplanes.  Each  of  the 
airplane's  differences  were  written  in  their  FCN.  Thus  the  result  was  in  a  form  which  only 
the  computer  could  interpret,  unlike  the  ship  image  summaries  prepared  by  the  ship 
analysis  program  presented  in  this  thesis. 

3 .  Identification  of  Spherical  Objects 

Computers  can  rapidly  process  large  amounts  of  data  but  have  difficulty 
identifying  objects  because  of  the  large  number  of  possibilities.  Humans  skilled  in 
photographic  interpretation  have  problems  processing  large  amounts  of  imagery  data,  but 
are  very  skilled  at  identifying  an  object  once  it  is  detected.  A  more  human-like  computer 
identifier  was  presented  by  Cox,  et  al.  [Ref.  3],  in  a  detector  for  all  spherical  objects 
contained  in  an  image. 

The  first  step  of  their  detection  methodology  is  parameter  extraction,  where  the 
position  of  the  light  source  and  image  intensity  parameters  are  gathered.  The  next  step  is 
feature  detection,  using  either  input  from  a  previous  program  or  by  conducting  a  simple 
pattern  search,  which  takes  about  60  seconds  for  a  512  by  512  image.  The  third  step  is 


segmentation,  which  uses  the  gradient  angle  transform  to  segment  an  image,  creating  a 
candidate  region  to  examine.  The  final  step  is  validation,  in  which  the  region  is  processed 
and  classified. 

Their  approach  limits  the  universe  of  objects  the  computer  must  identify  to  one 
type  while  using  several  different  detection  algorithms  to  ensure  detection.  The  approach 
of  this  thesis  also  limits  its  universe  significantly. 

4 .  Parallel  Processors  and  Image  Analysis 

Because  of  the  large  amounts  of  data  a  single  image  contains,  research  using 
parallel-processing  computers  is  also  being  conducted  [Ref.  6].  In  one  example 
implemented  using  a  Connection  Machine,  the  operator  identifies  an  object  for  the  machine 
to  find  and  gives  an  example.  The  program  studies  the  example  and  builds  a  database  of 
knowledge  about  the  object.  Using  this  knowledge,  it  then  examines  an  image,  finding  and 
highlighting  any  of  those  objects  present  in  the  image.  Using  the  Connection  Machine 
reduces  processing  time  from  30  minutes  to  under  one  second  for  pictures  ranging  in  size 
from  256  by  256  pixels  to  5000  by  5000  pixels  [Ref.  7]. 

This  approach  works  at  a  lower  level  of  analysis  than  the  approach  in  this  thesis. 
It  could  help  the  analyst  in  the  initial  processing  stage  to  highlight  areas  which  should  be 
examined  first. 

C.    RESEARCH  USING  AI  TECHNIQUES 
1 .  Knowledge-Based   Airport   Identification 

One  example  of  research  combining  conventional  and  knowledge-based  techniques 
is  SPAM,  system  for  airport  photographic  interpretation  using  MAPS  [Ref.  8].  Using 
conventional  image  processing  tools,  SPAM  first  detects  and  labels  regions.  Then,  using 
rule-based  control  and  recognition  methods,  it  groups  major  components  of  airports. 


This  technique  initially  suffered  from  the  same  limitation  mentioned  in  section  B-3, 
in  that  its  universe  of  knowledge  was  very  limited,  in  this  case  to  airports.  Later  research 
[Ref.  9]  focused  on  ways  to  improve  the  knowledge  acquisition  process  for  other 
universes,  resulting  in  a  set  of  interactive  tools.  These  applications  concentrate  on 
identifying  all  occurrences  of  a  structure,  whereas  this  thesis  concentrates  on  the  pieces  of 
an  object  within  a  structure. 

2 .  Rule-Based  Methods  for  Image  Analysis 

Recent  work  attempted  to  analyze,  segment,  and  interpret  images  of  printed  circuit 
boards  and  satellite  images  of  ice  flows  [Ref.  10].  The  methods  developed  used  a 
combination  of  procedural  and  rule-based  programming.  The  process  was  divided  into 
three  levels:  low,  mid,  and  high.  Low-level  processing  involved  image  operations  such  as 
thresholding,  computation  of  gradients  and  non-maxima  suppression,  etc. 

The  mid-level  processing  involved  grouping  of  similar  edges  and  extracting 
symbolic  entities.  The  mid-level  processing  extracted  the  initial  symbolic  information, 
including  the  geometric  properties  of  each  linked  segment,  the  forward  and  backward 
neighbor  codes,  and  the  codes  or  labels  of  edges  passing  through  pixels.  Further  mid-level 
processing  joined  small  segments  and  edges  together. 

The  high-level  processing  attempted  to  identify  the  extracted  entity,  deriving  a 
symbol  that  has  meaning  to  either  humans  or  another  follow-on  process.  This  step  used  a 
hypothesis  generation/reduction/verification  model.  The  goal  was  to  produce  a  symbolic 
description  of  the  input  image.  For  a  circuit  board,  this  would  be  a  description  of  the 
length,  location,  and  orientation  of  segments. 

The  program  works  towards  the  same  goal  as  ours  in  trying  to  take  a  low-level 
group  of  data  and  bring  the  description  of  the  object  to  a  symbolic  level  understood  by 
humans.    The  final  symbolic  outputs  of  the  programs  are  different.    Theirs  prints  the 


location  and  size  of  an  object,  while  ours  goes  a  step  further  and  identifies  what  type  of 
object  is  being  examined. 

3 .  Recognizing  Man-Made  Objects  in  Aerial  Images 

Another  system  for  analyzing  aerial  imagery  was  developed  entirely  with  AI 
languages,  specifically  Lisp  and  ART  (the  Automated  Reasoning  Tool)  [Ref.  11].  The 
approach  differed  from  [Ref.  10]  in  that  it  was  goal-driven:  the  system  identified  all 
requested  object  types  in  an  image. 

The  application  was  tested  with  airport  scenes  and  also  as  a  target-cueing  aid  for 
FLIR  (Forward -Looking  Infrared)  imagery.  For  low-level  processing,  the  system  could 
utilize  several  different  image-processing  techniques.  The  program  could  select  the  proper 
technique  based  on  parameters  such  as  the  application  area,  sensors  used,  resolution, 
weather  conditions,  and  quality  of  data.  The  low-level  processing  concentrated  on  edge 
detection,  region  segmentation,  and  other  image  enhancement  functions. 

Intermediate-level  processing  was  responsible  for  creating  and  filling  slots  for  each 
region.  These  slots  included  coordinates,  area,  variance,  connectivity,  compactness,  and 
other  descriptions. 

Knowledge  about  the  area  of  interest  improved  the  interpretation  process.  The 
system  included  information  on  basic  terrain  features  such  as  power  plants,  roads, 
railroads,  rivers,  forests,  etc.  Based  on  the  goal  given  the  system,  it  would  also  generate  a 
dynamic  model  for  the  goal  objects.  Then  the  high-level  processing  would  classify  regions 
into  goal  objects,  using  the  information  in  the  region  slots  while  attempting  to  make  the  best 
matches  possible. 

The  approach  provides  larger  flexibility  than  the  one  presented  in  this  thesis.  Also, 
instead  of  identification,  it  concentrates  on  detecting  and  highlighting  objects  for  the  human 
operator  to  examine. 


III.    METHODS  OF  SHIP  IDENTIFICATION 

For  people  to  develop  the  necessary  expertise  in  ship  identification  is  a  time- 
consuming  and  personnel-intensive  task.  There  are  hundreds  of  types  of  ships  to  be 
studied,  each  with  its  own  set  of  identifying  characteristics.  Learning  these  characteristics 
requires  days  of  study,  as  well  as  frequent  refresher  training  to  maintain  the  acquired  skills. 

Instruction  methods  for  image  interpretation  generally  follow  a  rule-based  approach. 
The  student  must  develop  his  working  knowledge  of  ships  and  then  memorize  key 
differences  between  types  of  ships.  The  first  step  is  to  teach  the  student  the  different 
features  of  a  ship.  These  features  usually  stand  out  from  a  silhouette  or  image  and  can  be 
peeled  from  the  ship  and  analyzed  independently.  During  this  analysis  the  feature's  identity 
is  further  refined  to  as  much  detail  as  possible.  This  is  shown  in  the  table  below. 


TABLE  3.1.     Example  Identifications  of  Ship  Structures 

Feature 

Intermediate  Identity 

Exact  Identity 

gun-turret 

1  5-in  gun  tube,  located 
forward,  square  turret 

Mark  54 

weapons-system 

square       launcher, 
forward,  several  tubes 

ASROC  Launcher 

radar 

large,  square  antenna, 
probably  air-search. 

Lockheed  SPS-40 

large  structure  aft 

box-like  shape,  clear 
area  beside,  several  antennas 
on  box 

Hangar  for  LAMPS 

Let  us  assume  the  student  is  learning  about  air-search  radar.  The  student  must  first 
learn  a  series  of  broad  rules  for  each  feature,  like  that  air-search  radar  usually  have  large, 
squared  antennas,  rotate  360  degrees,  and  are  located  high  in  the  ship  superstructure. 
Partial  classification  (such  as  recognition  that  the  feature  is  a  radar)  of  a  feature  can  allow 


the  ship  identification  process  to  still  proceed.  The  next  step  is  learning  to  distinguish 
between  different  air-search  radar.  For  example,  the  Lockheed  SPS-40  radar  antenna  has  a 
triangular  base,  a  square  antenna  with  tapered  ends,  and  a  feed  horn  that  hangs  over  the  top 
of  the  antenna. 

Once  a  particular  radar  has  been  accurately  identified,  the  student  must  learn  to 
associate  it  with  particular  ships.  For  example,  the  SPS-40  radar  is  found  on  several  US 
FFG  class  ships,  including  the  BRONSTEIN  and  GLOVER  classes.  The  student  then  tries 
to  learn  every  ship-class  known  to  use  the  SPS-40  radar.  Such  detailed  knowledge 
deteriorates  rapidly  without  frequent  use  or  refresher  training. 

Once  the  student  is  capable  of  identifying  the  different  features,  he  next  learns  the  rules 
for  identifying  particular  ships.  This  requires  learning  the  differences  that  the  same  features 
may  have  on  different  ships.  These  differences  include  the  way  different  features  are 
combined  as  well  as  where  on  the  ship  these  features  are  located.  The  learning  process  is 
again  rule-based,  as  in  the  GARCIA-class  FFG's  which  have  one  Mk-54  gun-turret 
forward  and  one  Mk-54  gun-turret  aft,  an  ASROC  launcher  behind  the  forward  turret,  a 
Lockheed  SPS-40  radar,  a  single  stack,  and  an  SPS-10  radar.  The  difference  between  the 
SAMPLE  and  BRADLEY,  both  GARCIA-class  FFG's,  is  the  BRADLEY  is  fitted  for  the 
LAMPS  helicopter  and  has  a  hanger  on  the  rear  deck.  The  rules  are  learned  for  all  classes 
of  ships,  and  the  skills  will  rapidly  deteriorate  without  constant  use. 


IV.     A  PICTURE-DESCRIPTOR  EXTRACTION  PROGRAM 

A.   OVERVIEW 

Our  feature  extraction  and  identification  process  involves  three  stages;  (1)  finding  the 
silhouette  boundary  (outline),  (2)  locating  the  bumps  on  the  silhouette,  and  (3)  extracting 
and  identifying  the  features  of  the  bumps.  Stage  1  is  implemented  in  C  and  stages  2  and  3 
are  written  in  MPROLOG.  The  program  flow  is  shown  in  Figure  4.1  below.  This  chapter 
is  broken  into  three  corresponding  sections. 


5hip_ 
Trace 
Program 

> 

f 

Locate 
Bumps 
Program 

> 

1 

Identify 

Bumps 

Program 

Figure  4.1.    Diagram  of  Program  Modules 
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Figure  4.2  contains  a  labeled  ship  silhouette,  demonstrating  some  bump  identifications 
our  program  attempts  to  make. 


radar      mast/support 
weapons-system  \^^  antenna 

J 


gun-tun 


uperstructure 


Figure  4.2.     Silhouette  with  Labeled  Bumps 


B  .   FINDING  THE  SILHOUETTE  BOUNDARY 

The  ship  silhouettes  are  digitized  using  an  Eikonix  digitizer  and  SUN  Workstation. 
The  digitizing  program  is  scanit,  written  by  David  S.  Hill  for  a  computer  graphics  class. 
The  ship  must  be  positioned  upright  and  level  with  the  bow  on  the  left  side.  The 
orientation  is  important  as  the  second  stage  exploits  the  identity  of  y-coordinates.  Using  an 
integration  value  of  about  12,000,  the  gitized  image  is  captured  into  a  file  as  a  1000  by 
1000  array  of  the  pixel  values.  This  file  is  transferred  to  the  VAX  11/785  computer  for 
boundary  processing. 

The  ship  silhouettes  used  were  taken  from  Jane's  All  the  World's  Fighting  Ships 
[Ref.  12].  Six  different  ships  of  the  frigate  class  were  used.  The  selection  criteria  for  the 
silhouettes  were  they  have  many  mutual  features  in  common,  with  differences  in  the  size 
and  location  of  the  features.  Each  digitized  silhouette's  file  is  processed  separately  by  the 
ship_trace  program,  with  the  resulting  files  of  line  segments  passed  on  to  the  next  stage. 


The  first  problem  to  be  solved  when  finding  the  boundary  of  the  digitized  image  is 
selecting  a  proper  threshold.  The  ship  silhouettes  were  solid  black  shapes  on  white  paper. 
Although  black  would  be  indicated  by  a  pixel  value  near  0  (0  hex)  and  white  would  be  near 
255  (ff  hex),  setting  the  threshold  to  catch  the  change  from  black  to  white  requires  some 
adjustment  for  each  picture.  This  was  necessary  to  ensure  that  any  noise  captured  with  the 
image  was  tuned  out.  Also,  a  threshold  too  close  to  black  would  tend  to  fill  in  many  of  the 
fine  details  of  the  superstructure,  while  a  threshold  too  close  to  white  would  allow  breaks 
in  antennas  or  cracks  through  the  ship  to  appear  in  the  image.  To  help  in  selecting  the 
proper  threshold  value,  a  small  program  was  written  which  printed  out  a  portion  of  the 
silhouette  with  asterisks  for  the  dark  areas  and  spaces  for  the  light  ones.  By  running  this 
program  using  different  values  and  viewing  the  results,  the  threshold  value  could  be  fine- 
tuned  to  capture  the  right  amount  of  detail. 

Once  the  threshold  value  is  selected,  the  boundary-finding  program  can  be  run.  This 
program  was  designed  to  walk  along  the  boundary1,  of  the  silhouette,  recording  (x,y) 
coordinates  of  points  whenever  it  changed  direction.  The  program  uses  the  four  principal 
directions,  up,  down,  left,  and  right.  Its  goal  was  to  start  at  the  left  end  and  keep  moving 
toward  the  right  end  of  the  ship.  The  approach  is  similar  to  Pavlidis'  contour  tracing 
algorithm  [Ref  13:p.  143].  The  program  in  effect  traces  the  outline  of  the  ship,  creating  a 
file  of  turn  points.  Each  pair  of  points  denotes  a  line  segment.  A  long  straight  antenna 
would  be  converted  into  two  or  three  points,  while  a  gun  turret  might  create  fifty  or  more 
points. 


^e  boundary  is  located  along  the  black  cells  delineating  the  top  of  the  ship.  In  effect, 
the  program  would  walk  the  deck  of  the  ship,  recording  locations  where  it  changed 
direction.  The  locations  were  expressed  in  the  coordinates  of  the  black  cells. 
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Originally,  we  intended  to  trace  the  oudines  using  eight  directions,  the  four  explained 
above  as  well  as  the  diagonals.  But  diagonal  lines  can  also  be  represented  by  a  pattern  of 
horizontal  and  vertical  moves.  Although  this  creates  many  more  turn  points,  finding  and 
recording  diagonals  requires  looking  at  a  larger  group  of  points  during  each  iteration  than 
our  algorithm  was  designed  for. 

Turn  points  found  are  written  to  a  data  file.  This  data  file  is  then  converted  to  an 
acceptable  MPROLOG  format:  a  list  with  each  point  in  the  list  being  represented  as 
cp(<x-coordinate>,<y-coordinate>). 

C.   LOCATING  THE  BUMPS 

The  extraction  algorithm  was  designed  to  look  for  matching  pairs  of  up-turns  and 
return-to-horizontals.  These  are  the  definers  of  a  single  "bump",  the  point  where  it  starts 
going  up  and  the  point  where  it  returns  to  the  horizontal  at  a  height  equal  to  that  of  the  up- 
turn. By  detecting  those  two  points,  and  capturing  all  the  line  segments  between  them,  a 
bump  can  be  extracted.  A  modified  algorithm  is  used  to  find  bumps  which  start  and  end  at 
different  heights.  The  program  repeatedly  examines  the  point  where  it  currently  is,  the 
previous  point,  and  the  next  point.  From  these,  it  determines  if  the  middle  point  was 
located  at  an  up-tum.  If  so,  it  saves  the  current  location  coordinates  in  a  list  of  up-turns. 
This  is  repeated  for  every  successive  triple  of  points. 

In  order  to  capture  only  true  comers,  an  up-turn  was  required  to  have  at  least  a  three- 
pixel  vertical  movement.  This  helps  ignore  parts  of  the  silhouette  which  were  jagged  due  to 
noise.  The  retum-to-horizontal  which  matched  these  up-turns  was  only  required  to  have  a 
one-pixel  movement.  By  having  the  up-tum  constraint  select  only  well-defined  corners, 
the  retum-to-horizontal  constraint  could  be  looser  to  make  it  easier  to  find  the  matching 
comer  of  the  bump.  Figure  4.3  shows  examples. 


3  pixel  up_turn 


\ 


*■  ■*  * 
*  •*  ■* 


Matching  return_to_horizontal 
Noise 


*  *  *       /  Noise 

»fl/       /|\ 


*  *  *■  •*  *  *[*  *  * 


Figure  4.3.    Example  of  Normal  Bump  and  Image  Noise 

Not  all  bumps  on  a  ship  can  be  caught  this  way,  so  two  variations  of  the  algorithm  are 
required.  One  operates  as  an  inverse  of  the  first.  If  the  return-to-horizontal  has  a  three- 
pixel  or  more  vertical  movement  preceding  it,  the  program  turns  around  and  retraces  the 
bumps  until  it  finds  a  one-pixel  or  two-pixel  vertical  movement  where  the  y  coordinates  are 
equal  to  the  retum-to-horizontal's.  An  example  is  shown  in  Figure  4.4. 


*  *          3  pixel 
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return_to_horizontal 

Created  up_turn 
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Figure  4.4.     Example  of  Created  Up-turn 
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The  final  variation  is  designed  to  detect  and  create  a  bump  when  a  vertical  movement 
passes  below  a  previously  labeled  up-turn.  That  terminates  all  unmatched  bumps  whose 
up-turns  were  in  the  vertical-movement  height  range.  An  example  is  shown  in  Figure  4.5. 


3  pixel  up_turn 
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********************* 


Figure  4.5.    Example  of  Created  ReturntoHorizontal 

The  file  containing  matching  pairs  of  up-turns  and  return-to-horizontals  is  kept  in 
MPROLOG  format.  When  the  algorithm  reaches  the  end  of  the  ship,  it  closes  off  any 
bumps  still  in  the  up-turn  list  and  closes  the  file. 

D.    BUMP  FEATURE  IDENTIFICATION 

1 .  Removing  the  Shape  from  the  Silhouette 

The  shape  analysis  program  is  written  in  MPROLOG.  It  iterates  through  a  list  of 
bumps,  finding  their  size,  length-to-width  ratio,  curviness,  and  location.  These  descriptors 
are  then  assigned  categories  and  rules  for  classification  of  bumps  from  descriptor  categories 
are  applied. 

To  do  this,  the  program  traces  the  outline  of  the  bump  using  the  list  of  turn  points. 
If  the  bump  has  other  bumps  within  it,  these  bumps  are  stripped  off  so  only  the  basic  shape 
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remains.    If  a  bump  within  the  bump  is  removed,  the  program  creates  a  straight  line 
segment  where  the  bump  was.  This  is  shown  in  Figure  4.6. 


peeled  bump 

*  * 

*  *  *  *  * 
*  *  x  *  x  * 

■*  ■*  ■*  ■* 

J 

x  *  *  „  „  „  n  \m 
*•  *  *■  ■*  *•  x  *  1 

| inserted  line  segment 

1*  *•  *•  •*  *  * 

*■  *  *•  ■* 
*  *■  *  * 
■*  *-  ■*  •* 

~ simplified  bump 

Figure  4.6.     Example  of  Peeling  Bump  and  Creating  Line  Segment 

Next,  each  bump  is  analyzed  as  to  curviness,  horizontal  location,  length-to-width  ratio, 
orientation,  and  size.  The  number  of  turn  points  in  the  bump  list  roughly  indicates  the 
curviness  of  the  shape  (See  Figure  4.7).  The  location  of  the  shape  is  determined  as  bow, 
forward,  mid_ship,  aft,  or  stem;  this  is  useful  as  gun  systems  tend  to  be  forward  or  aft  on 
ships,  while  radar  systems  are  usually  found  in  the  middle  third  of  the  ship,  around  the 
superstructure.  The  length-to-width  ratio  of  the  object  is  found  from  the  minimum  and 
maximum  x  and  y  values,  which  create  an  imaginary  box  around  the  object,  as  exemplified 
in  Figure  4.8.  Particular  length-to-width  ratios  can  be  described  as  the  shapes  square, 
pole,  rectangle,  etc.  Box  orientation  (flat_rectangle,  tall_rectangle)  can  also  be  described. 
The  box  size  in  relation  to  the  ship  is  also  used. 
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Actual  Curve 
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Figure  4.7.     Example  of  Curve  Approximation 
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Figure  4.8.     Example  of  Box  Surrounding  a  Shape 
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2.  Shape  Analysis  and  Identification 

After  all  descriptors  of  a  bump  have  been  found,  shape  rules  (for  antenna,  radar, 
gun-turret/weapons-system,  weapons-system,mast/support,or  superstructure)  are  applied. 
For  example,  the  descriptors  [pole,  tall_rectangle,  medium,  mid_ships]  indicate  an  antenna 
of  some  type,  and  "antenna"  is  put  into  the  identification  list 

Sometimes  a  bump  will  match  more  than  one  rule.  For  example,  a  shape  may  have 
attributes  of  a  radar  as  well  as  those  of  a  superstructure.  Then  both  rule  descriptions  are 
added  to  the  list  and  the  word  "or"  is  inserted  between  them. 

Sometimes  a  bump  will  not  match  any  of  the  rules.  This  may  be  due  to  noise 
introduced  when  the  image  was  digitized,  or  the  bump  might  be  an  uncommon  shape  for 
these  types  of  ships.  In  this  case,  the  word  "unknown"  is  inserted  in  the  identity  list. 
Since  the  location  of  the  shape  and  its  description  list  is  also  kept  for  future  use,  a  more 
intelligent  program  might  later  be  able  to  identify  the  shape. 

E.    FINAL  OUTPUT  OF  THE  PROGRAM 

When  the  programs  have  completed  all  processing,  the  final  file  consists  of  a  listing  of 
all  the  extracted  shapes,  their  descriptors,  and  their  identity,  if  known.  The  shapes  are 
identified  by  their  starting  and  ending  coordinates,  their  descriptor  list,  and  their 
identification  list.  This  file  is  then  available  for  processing  by  a  future  program  which 
could  match  a  ship  class  to  the  descriptors  our  program  has  created. 


V.  RESULTS  OF  PICTURE-DESCRIPTOR  EXTRACTION 
PROGRAM 

A.   FINDING  THE  SILHOUETTE  BOUNDARY 

The  performance  of  the  first-stage  program  is  shown  in  Table  5.1  below.  The 
measurements  indicate  the  size  of  the  MPROLOG  file  produced,  the  number  of  turn  points 
in  the  file,  and  the  threshold  setting. 


TABLE  5.1.     Outline  Extraction  Program  Results 

Ship  Name 

MPROLOG  File 
Size 

Number  of 
Points 

Threshold 

GARCIA  (LAMPS) 

12.5  Kb 

962 

60 

GLOVER 

17.9  Kb 

1,380 

50 

BROOKE 

12.4  Kb 

952 

50 

BRONSTEIN 

8.1  Kb 

634 

50 

KNOX 

10.4  Kb 

798 

55 

GARCIA 

11.9  Kb 

913 

50 

As  the  table  shows,  the  resulting  MPROLOG  file  is  not  very  large,  since  each  of  the 
original  digitized  image  files  was  990  Kb  in  size.  The  only  problem  encountered  involved 
the  silhouette  of  the  BRONSTEIN:  its  digitized  image  had  more  random  noise  than  the 
others,  and  the  program  required  several  runs  to  find  the  bow  of  the  ship,  the  starting 
point. 
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B  .  LOCATING  THE  BUMPS 

The  second  stage  of  finding  and  marking  the  location  of  bumps  on  the  silhouettes  was 
the  slowest  stage.  The  program  could  process  any  ship  except  the  GLOVER  in  one  run 
(the  number  of  turn  points  in  the  GLOVER's  data  file  had  to  be  split  in  half  to  allow  the 
program  to  run,  for  otherwise  the  program  crashed  after  reading  in  the  data  points  and  the 
required  program  modules).  A  listing  of  the  CPU  time  used  is  shown  in  Table  5.2. 


TABLE  5.2.    Execution  Time  for  Bump  Locating  Program 

Ship  Name 

CPU  Time  (minutes:seconds) 

GARCIA  (LAMPS) 

31:07 

GLOVER 

39:11 

BROOKE 

30:36 

BRONSTEIN 

13:17 

KNOX 

21:37 

GARCIA 

27:47 

Table  5.3  summarizes  the  number  of  bumps  found. 
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TABLE  5.3.     Results  of  the  Bump  Locating  Program 

Ship  Name 

Image  Length     (Pixels) 

Number  of  Bumps 

GARCIA  (LAMPS) 

735 

37 

GLOVER 

933 

48 

BROOKE 

769 

44 

BRONSTEIN 

700 

34 

KNOX 

757 

37 

GARCIA 

758 

26 

To  test  the  program  with  a  different  size,  the  GLOVER  silhouette  was  digitized  so  that 
it  completely  filled  the  imaging  window,  adding  more  detail  to  the  image  than  the  others, 
and  this  silhouette  generated  the  largest  number  of  bumps.  The  smallest  number  of  bumps 
was  generated  by  the  GARCIA-class  FFG  without  a  LAMPS  silhouette  because  the  rear  of 
the  ship  has  few. 

C.   BUMP  FEATURE  IDENTIFICATION 

Table  5.4  summarizes  the  results  of  the  third  program  in  the  identification  process.  A 
successful  identification  is  when  a  bump  is  identified  either  as  itself  or  as  a  part  of  an  or'd 
list.  If  the  bump  is  classified  as  unknown,  it  is  not  counted  as  a  successful  identification. 
The  use  of  unknown  as  a  descriptor  prevents  the  program  from  a  wrong  identification,  i.e. 
identifying  a  radar  as  a  gun-turret 
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TABLE  5.4.     Statistics  of  Bump  Identification  Program 

Ship  Name 

Number 
Identified 

Number 
Unknown 

Number 
Not  Done 

Percent 
Successful 

GARCIA 

(LAMPS) 

32 

4 

1 

86% 

GLOVER 

40 

8 

0 

83% 

BROOKE 

38 

6 

0 

86% 

BRONSTEIN 

26 

8 

0 

76% 

KNOX 

27 

10 

0 

73% 

GARCIA 

14 

6 

6 

53% 

TOTALS 

177 

42 

7 

78% 

As  the  table  shows,  the  identification  program  is  successful  in  identifying  a  bump 
about  three  of  four  tries.  The  number  of  bumps  a  ship  has  to  extract  does  not  seem  to 
affect  the  program's  success  rate.  A  bump  falling  in  the  unknown  category  is  not 
necessarily  a  failure;  it  may  have  an  unusual  shape,  there  may  be  noise  in  the  image  that  has 
distorted  its  shape,  or  the  detail  from  the  silhouette  may  not  provide  enough  information  to 
make  a  definitive  identification.  Some  features  are  more  helpful  in  the  identification 
process  than  others.  For  frigate-class  ships,  the  key  features  are  the  number,  location,  and 
type  of  weapons  systems,  the  number  of  radar  and  antennas,  and  the  presence  or  absence 
of  certain  superstructures  like  helicopter  hangars. 

To  illustrate  bump  features,  the  differences  between  a  GARCIA-class  FFG  and  the 
GARCIA  FFG  with  LAMPS  can  be  seen  in  Figure  5.1. 
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GARCIA-class  FFG 


GARCIA  FFG  vith  LAMPS 


Figure  5.1.     Silhouettes  of  GARCIA-class  FFGs 

The  ships  have  the  same  basic  structure  and  similar  weapons  systems  and  radar.  The 
key  differences  are  the  number  of  antennas  and  the  presence  of  the  helicopter  hangar  on  the 
rear  deck  of  the  FFG  with  LAMPS.  Table  5.5  describes  a  GARCIA-class  FFG  without  the 
LAMPS  helicopter  modification.  Comparing  it  and  Table  5.6  illustrates  the  differences  our 
program  detects  between  two  similar  ship  silhouettes. 
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TABLE  5.5.     Features  for  GARCIA-class  FFG  (no  LAMPS) 

Feature  on  the  Ship 

Feature  Identified  by  Program 

gun-turret 

weapons-system 

weapons-system 

weapons-system 

radar 

radar  or  superstructure 

radar 

radar  or  superstructure 

radar 

(program  cannot  process) 

weapons-system 

weapons-system 

antenna 

antenna 
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TABLE  5.6.    Features  for  GARCIA  FFG  with  LAMPS 

Feature  on  the  Ship 

Feature  Identified  by  Program 

gun-turret 

gun-turret/weapons-system 

weapons- system 

weapons- system 

antenna 

antenna 

antenna 

antenna 

radar  and  mast 

mast/support 

antenna 

antenna 

mast/s_up2ort 

mast/support  or  superstructure 

antenna 

antenna 

antenna 

antenna 

gun-turret 

gun-turret/weapons-system 

antenna 

antenna 

antenna 

antenna 

hangar  superstructure 

superstructure 

D.   PROGRAM  LIMITATIONS  DUE  TO  HARDWARE  CONSTRAINTS 

A  ship's  bumps  could  be  identified  in  two  to  five  runs  of  the  program.  The 
requirement  for  multiple  runs  was  due  to  hardware  and  software  limitations  of  the  ISI 
Workstations.  The  MPROLOG  statement  table  does  not  use  any  garbage-collection 
techniques  to  free  up  memory  after  it  is  no  longer  needed.  Because  of  the  large  amount  of 
recursive  list  processing  used  in  the  extraction  process,  the  statement  table  fills  up  quickly 
and  the  program  halts.  After  moving  the  coordinates  of  bumps  that  have  been  processed  to 
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the  end  of  the  list  of  bumps,  the  program  can  be  restarted.  It  will  continue  to  process  the 
remaining  bumps,  adding  their  descriptions  as  they  are  developed. 

A  similar  limitation  experienced  by  the  feature  extraction  and  identification  program 
was  its  inability  to  process  seven  bumps  when  analyzing  the  two  GARCIA-class  frigates. 
The  error  occurs  as  the  program  attempts  to  strip  off  several  bumps  mounted  on  top  of  the 
ships'  central  superstructure.  The  number  of  points  to  be  extracted  for  the  bump  and  the 
number  of  points  being  processed  to  peel  off  a  bump  are  too  many  for  the  system  stack. 
However,  if  the  feature  causing  the  problem  is  removed  from  the  list  and  the  program  is 
restarted2  ,  the  program  operates  normally.  The  effects  of  this  inability  to  process  six 
bumps  is  shown  in  the  shorter  length  of  Table  5.5  and  the  overall  lower  identification 
success  rate  for  the  GARCIA-class  FFG  without  LAMPS.  The  problem  can  be  solved 
either  by  adding  more  memory  or  changing  the  software. 


2  In  this  case,  the  ordering  of  the  list  of  bumps  is  changed,  either  moving  the  problem 
bumps  behind  the  end  of  list  flag  or  by  removing  them  from  the  list.  The  program  is  then 
restarted  from  the  beginning,  unlike  the  restarts  mentioned  earlier  where  the  program  picks 
up  where  it  left  off. 
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VI.     CONCLUSION 

A.   SUCCESS  OF  FEATURE  EXTRACTION  AND  IDENTIFICATION 

We  have  demonstrated  that  the  key  features  from  a  ship  silhouette  can  be  extracted, 
analyzed,  and  identified  to  produce  a  list  describing  and  naming  the  features.  The  results 
we  have  seen  are  similar  to  other  research  using  knowledge-based  and  artificial-intelligence 
techniques.  The  average  success  rate  for  six  frigate-class  ships  was  78%,  with  values 
ranging  from  53%  to  86%. 

These  programs  provide  a  versatile  front-end  method  for  automated  ship  recognition 
and  identification.  The  description  and  identification  information  contained  in  the  final 
listing  could  be  used  by  other  programs  to  identify  the  type  of  ship  or  possibly  the  name  of 
the  ship.  Such  a  program,  combined  with  an  intelligent  tutor  program,  would  relieve 
instructors  of  the  administrative  and  repetitive  tasks  of  updating  silhouette  libraries. 

The  MPROLOG  programs  suffer  from  hardware  and  software  limitations.  The 
amount  of  information  contained  in  a  silhouette  is  not  that  large,  but  the  recursive  analysis 
used  by  the  programs  quickly  use  up  available  memory  and  stack  space.  Many  other 
Prolog  dialects  should  not  have  these  limitations. 

B  .   AREAS  FOR  FUTURE  WORK 

The  program  could  be  extended  to  handle  all  ships,  regardless  of  nationality  or  type. 
The  program  could  also  be  improved  by  adding  the  ability  to  detect  and  analyze  objects 
with  diagonal  appendages.  Our  program  converts  diagonals  into  horizontal  and  vertical 
steps.  Extending  it  to  interpret  diagonals  would  allow  the  program  to  better  pick  apart 
pieces  of  the  superstructure  and  identify  more  of  the  antenna  and  mast  detail. 
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Adding  the  ability  to  reason  from  known  bumps  would  allow  the  program  to  go  back 
and  reevaluate  bumps  previously  identified  as  unknown  and  make  an  identification. 
Identification  is  possible  by  knowing  what  ship  classes  partially  match  the  current 
silhouette  and  looking  for  unknown  bumps  that  could  match  known  ship  structures  in  the 
ship  classes  being  examined. 

These  programs  could  be  put  together  with  an  intelligent  tutoring  program.  The 
combined  system  would  be  a  useful  tool  for  developing  and  maintaining-silhouette 
identification  skills. 
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APPENDIX  A  -  C  PROGRAM  SOURCE  CODE 


SHIP_TRACE.C 

This  program  finds  the  bow  and  stern  of  the  digitized  image  and  creates  a  binary  file 
line. coords  is  then  used  by  program  format. c  to  create  a  MPROLOG  data  file  of  the 
turn  points. 


#include  <stdio.h> 

main() 
/ 

l 

char 

PTS  [250][1000]; 

char 

row  data[1000]; 

int 

num_pixels[1000]; 

int 

datafile,pfile,pfile  1 ; 

int 

XCOORD,YCOORD; 

char 

code; 

int 

X,Y; 

int 

row,i; 

int 

fd[l]; 

int 

info[2j; 

short 

TYPE,ROWS  .COLUMNS; 

int 

max_rows,  max_columns,  max_PTS; 

int 

xstartpoint,  ystartpoint; 

int 

xendpoint,  yendpoint; 

int 

threshold;                                /*  adjust  to  picture  quality  */ 

/*    open  datafile    */ 
printf ("starting  program  \n"); 

if  ((datafile  =  open(7scratch/bizer/pic.brooke'',0))  <  0  ) 
{ 
printf("can't  find  it\n"); 
exit(); 
} 

/*    open  output  file  */ 
if  ((pfile  =  creat("line.coords",0744))  <  0  ) 
{ 

printf("can't  open  output  line.coords\n"); 

exit(); 


/*     open  output  file  */ 
if  ((pfilel  =  creat("line.info",0744))  <  0  ) 
{ 

printf  ("can't  open  output  line,  info  \n"); 

exitQ; 
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} 

/*  pull  off  header  information  */ 

read(datafile,&TYPE,2); 

read(datafile,&ROWS,2); 

read(datafile,&COLUMNS,2); 

printf("type  %d  max_rows  %d  max.columns  %d\n",TYPE,ROWS,COLUMNS); 

max_rows= 1 000; 

max_columns  =  1000; 

/♦clean  off  empty  space  above*/ 
for  (row  =  0;  row  <  280;  row++) 

read(datafile,  &row_data[0],  max_columns); 
printf  ("cleaned  off  blank  space  \n"); 

/*    get  a  row  of  information    */ 
for  (row  =  280;  row  <  510;  row++) 

read(datafile,  &PTS[row-280][0],  max_columns); 

printf("  Array  is  initializedAn"); 

/******* ******************************************************* ********/ 
/*    Initialize  tracing  program  */ 

xstartpoint  =  1000; 
xendpoint  =  0; 

for  (row  =  0;  row  <  220;  row  ++) 
for  (i  =  20;  i  <  max_columns;  i++) 
{ 

if  (PTS[row][i]  <  55  &&  i  <  xstartpoint) 
{ 
xstartpoint  =  i; 
ystartpoint  =  row; 


if  (PTS[row][i]  <  55  &&  i  >  xendpoint) 
{ 

xendpoint  =  i; 

yendpoint  =  row; 


XCOORD  =  xstartpoint; 
YCOORD  =  ystartpoint; 
X  =  XCOORD; 
Y  =  YCOORD +  280; 

ystartpoint  =  ystartpoint  +  280; 
yendpoint  =  yendpoint  +  280; 

printf("XStartpoint  =  %d  XEndpoint  =  %d\n",xstartpoint,xendpoint); 
printf("YStartpoint  =  %d  YEndpoint  =  %d\n",ystartpoint,yendpoint); 
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write(pfile  1  ,&y  startpoint,4); 
write(pfile  1  ,&xstartpoint,4); 
write(pfile  1  ,&yendpoint,4); 
write(pfile  1  ,&xendpoint,4); 
printf("made  it  past  writing  to  line.info\n"); 

code  =  'R'; 

prev_code  =  Z'; 

threshold  =  55; 

while  (XCOORD  <  xendpoint) 

{ 

while  (code  ==  'R') 

{ 
if  (PTS[YCOORD-l][XCOORD]  <  threshold) 

( 

printf("entering  with  code  =  R\n"); 

write(pfile,&Y,4); 

write(pfile,&X,4); 

YCOORD=  YCOORD-1; 

Y  =  Y+1; 
code  =  'U'; 
} 

else  if  (PTS[YCOORD][XCOORD+l]  <  threshold) 

{ 

X  =  X+  1; 

XCOORD =  XCOORD +  1; 

} 
else  if  (PTS[YCOORD+l]  [XCOORD]  <  threshold) 

{ 

write(pfile,&Y,4); 

write(pfile,&X,4); 

YCOORD  =  YCOORD  +  1; 

Y  =  Y-1; 
code  =  'D'; 
) 

else  if  (PTS[  YCOORD]  [XCOORD- 1]  <  threshold) 
{ 
printf( "entering  with  code  =  R\n"); 
write(pfile,&Y,4); 
write(pfile,&X,4); 
X  =  X-  1; 

XCOORD  =  XCOORD -1; 
code  =  'L'; 


while  (code  ==  'D') 
{ 

if  (PTS[YCOORD][XCOORD+l]  <  threshold) 
{ 
write(pfile,&Y,4); 


write(pfile,&X,4); 
XCOORD  =  XCOORD  +  1; 
X  =  X+1; 
code  =  'R'; 
} 
else  if  (PTS[YCOORD+l][XCOORD]  <  threshold) 

{ 

YCOORD  =  YCOORD  +  1; 

Y-Y-l; 

} 
else  if  (PTS[YCOORD][XCOORD-l]  <  threshold) 
{ 

write(pfile,&Y,4); 
write(pfile,&X,4); 

x  =  x  -  r 

XCOORD  =  XCOORD -1; 
code  =  'L'; 
} 
else  if  (PTS[YCOORD-l]  [XCOORD]  <  threshold) 
{ 

write(pfile,&Y,4); 
write(pfile,&X,4); 
YCOORD  =  YCOORD  -  1; 
Y  =  Y+1; 
code  =  'IT; 


while  (code  ==  'L') 
{ 
if  (PTS[YCOORD+l] [XCOORD]  <  threshold) 

{ 

write(pfile,&Y,4); 

write(pfile,&X,4); 

YCOORD  =  YCOORD  +  1; 

Y  =  Y-1; 
code  =  'D'; 
} 

else  if  (PTS[YCOORD][XCOORD-l]  <  threshold) 

{ 

X  =  X-  1; 

XCOORD  =  XCOORD- 1; 

} 
else  if  (PTS[YCOORD-l][XCOORD]  <  threshold) 

{ 

write(pfile,&Y,4); 

write(pfile,&X,4); 

YCOORD  =  YCOORD- 1; 

Y  =  Y+1; 
code  =  'U'; 
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else  if  (PTS[YCOORD][XCOORD+l]  <  threshold) 

{ 

write(pfile,&Y,4); 

write(pfile,&X,4); 

X  =  X+1; 

XCOORD  =  XCOORD  +  1; 

code  =  'R'; 


while  (code  ==  V) 
{ 
if  (PTS[YCOORD]  [XCOORD- 1]  <  threshold) 

{ 

write(pfile,&Y,4); 
write(pfile,&X,4); 
X  =  X  -  1* 

XCOORD  =  XCOORD -1; 
code  =  'L'; 
} 
else  if  (PTS[YCOORD-l]  [XCOORD]  <  threshold) 
{ 
YCOORD=  YCOORD-1; 

Y  =  Y+1; 
} 

else  if  (PTS[  YCOORD]  [XCOORD+1]  <  threshold) 

{ 

write(pfile,&Y,4); 

write(pfile,&X,4); 

XCOORD  =  XCOORD  +  1; 

X  =  X+  1; 

code  =  'R'; 

} 
else  if  (PTS[YCOORD+l][XCOORD]  <  threshold) 

{ 

write(pfile,&Y,4); 

write(pfile,&X,4); 

YCOORD  =  YCOORD  +  1; 

Y  =  Y-1; 
code  =  'D'; 


}  /*    end  while  X  <  endpt  */ 

close(datafile); 

close(pfile); 

close(pfilel); 

}/*    end  pgm    */ 
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FORMAT.C 

This  program  reads  line.coords,  converts  the  turn  points  to  MPROLOG  format  and 
writes  the  points  to  file  coords. pro. 

#include  <stdio.h> 
#include  <sys/file.h> 


mainO 
1 

FILE 

*datafile,  *datafilel,  *pfile; 

int 

row; 

int 

i; 

int 

xcoord,  ycoord; 

int 

max_rows; 

int 

testvalue; 

int 

first; 

/*  open  datafile    */ 

if  ((datafile  =  fopen("line.info",  "r"))  ==  NULL)  { 

printf("can't  find  it  infoW); 

exit(); 
} 
printf("opened  line.infoW); 

if  ((datafile  1  =  fopen("line.coords",  "r"))  ==  NULL)  | 

printf("can't  find  it  coords  \n"); 

exit(); 

} 
printf( "opened  line.coordsW); 

/**  open  output  file     **/ 

if  ((pfile  =  fopen("coords.pro","w"))  ==  NULL  )  { 

printf("can't  open  output  \n");  exit(); 

} 
printf( "opened  output  fileW); 

fprintf(pfile,"module  coords.W); 
fprmtf(pfile,7*$eject*^i  "); 
fprintf(pfil'  "bodyAn"); 
fprintf(pfik,"\n  "); 

/*  pull  off  starting  and  ending  coordinates  */ 
/*  coords  are  read  as  Y,X  and  written  as  X,Y  */ 

printf("before  fscan  \n"); 
fread(&ycoord,  sizeof(int),l, datafile); 
fread(&xcoord,  sizeof(int),  1  ,datafile); 
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fprintf(pfile,"start_pt([%d,%d]).\n",xcoord,ycoord); 

fread(&y  coord,  sizeof(int),l,datafile); 

£read(&xcoord,  sizeof(int),l,datafile); 

fprintf(pfile,"end_pt([%d,%d]).\n",xcoord,ycoord); 
fprintf  (pfile,"\npoints([ "); 
printf(" starting  to  process  rows\n"); 

I*  read  a  row  of  data  and  write  back  out  in  MPROLOG  format  */ 
for  (first=l;first  <  1999;first++)  { 

if  (fread(&ycoord,  sizeof(int),l,datafilel)  =  0) 
break; 

if  (fread(&xcoord,  sizeof(int),l,datafilel)  —  0) 
break; 

if  (first  !=  1) 

fprintf  (pfile,",\n"); 

fprintf(pfile,  "cp(%d,%d)",  xcoord,  ycoord); 
} 

fprintf(pfile;,]).\n"); 

fprintf(pfile,"endmod /*  coords.pro  */Nn"); 
fclose(datafile); 
fclose(datafilel); 
fclose(pfile); 
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APPENDIX  B  -  MPROLOG  PROGRAM  SOURCE  CODE 


PROCESS.PRO 

This  program  examines  the  turn  points  from  coords.pro  and  locates  all  bumps  on  the 
silhouette. 

module  process. 

/*$eject*/ 

body. 

import(add_statement/l,  del_statement/l). 

dynamic  (coordjist  /  1 ). 
dynamic  (cp  /  2). 
dynamic  (pts/  1). 
dynamic  (index_counter  /  1). 
dynamic  (up_turn  /  1). 
dynamic  (bump  /  4). 

totalJength(L)  :- 

start_pt([X,Y]),  end_pt([Xl,Yl]),  L  is  Xl-X  . 

set_up(L) :-  points(L), 
asserta(pts(L)), 
stars, 

write('START_POINT  -->,),start_pt([X,Y]),write(X),write('  '),write(Y), 
write(,END_POINT-->'),end_pt([Xl,Yl]),write(Xl),write('  '),write(Yl), 
nl,  stars, 

open(l, "bumps. pro"), 
tell(l, "bumps. pro"), 
write('module  bumps. '),nl, 
write(7*$eject*/,),nl, 
write('body.'),nl,nl, 
write('bumps(['), 
told(l). 

go :-    state(system_time,START), 
add_statement(up_turn([])), 
stars, 
set_up(L), 
run(X,Y), 
display_stats, 

state(system_time,FTNISH), 
REALJITME  is  FINISH  -  START, 
writeC Actual  Time  =  "),  write(R£AL_TIME), 
write("  shown  as  HHMMSS"). 
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run(X,Y)  :-  poims(L),indexl(L,0,PP,P,PN). 

get_pts(PP,P,PN\XP,YP,X.Y,XN,YN), 

walk(XP,YP^C,Y,XN,YN).add_statement(index_counter(0)), 

iterate(runl), 

close_bump_file. 

runl  :-  points(L),index_counter(K),K2  is  K+l,         f*  Iterates  this  pred  */ 
indexl(L,K2,PP.P,PN).  /*    until  runs  out    */ 

get_pts(PP,P,PN,XP,YP,X,Y,XN,YN),      /*    of  coordinates    */ 
walk(XP,YP,X,Y,XN,YN),  /*    to  process.       */ 

del_statement(index_counter(K)), 
add_statement(index_counter(K2)), 
display_stats,!. 

done  :-  points(L),length(L,N),index_counter(K),  K>=N-3. 

display_stats  :-  nl, stars. 

state(main_stack,  [U,C] ), 

write("main_stack  used  =  "),write(U),nl. 

state(statement_table,[Ul,Cl]), 

u-rite("statement_table  used  ="), 

write(Ul),write("  "),write(Cl),nl, 

state(cpu_time,TF), 

stars, writeC'cpu  time  =  "),write(TF),nl, 

stars. 

get_pts(PP,P,PN,XP,YP,X,Y,XN,YN)  :-  /*  Get  three  points  to  process  */ 

PP  =..  [cp,XP,YP], 
P  =..  [cp.X,Y]. 
PN  =..  [cp.XN.YN]. 


walk(XP,YP,X.Y,XN,YN)  :- 

change_up(XP.YP,X.Y,XN,YN). 
walk(XP,YP,X,Y,XN,YN)  :- 

below_upturn(XP,YP,X,Y,XN,Y.\). 
walk(XP,YP,X,Y,XN,YN)  :- 

rtn_horizontal(XP,YP,X.Y,XN,YN). 
walk(XP,YP,X,Y,XN,YN)  :- 

create_upturn(XP,YP,X,Y,XN,YNj. 
walk(XP,YP,X,Y,XN,YN)  :-  !. 

change_up(XP,YP,X,Y,XN,YX)  :-  XP<X,  /*  record  an  up_tum. 

L  is  YN  -  Y, 
L>2, 
up_tum(UL),append([(X,Y)],UL,ULl), 

add_statement(up_turn(ULl )), 
del_statement(up_tum(UL)),!. 
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rtn_horizontal(XP,YP,X,Y,XN,YN)  :- 
YP>Y 
Y=YN, 
X<XN, 
up_tum(UL), 
member((XU,Y),UL), 
delete((XU,Y),UL,ULl), 
del_statement(up_turn(UL)), 
add_statement(up_turn(UL  1 )), 
stars, 

write_to_screen(XU,Y,X,Y), 
write_to_file(XU,Y,X,Y),nl,!. 


/*  bump  rtn_to_horiz  at  same    */ 
/*   level  it  started,  so      */ 
/*    write  coordinates.        */ 


below_upturn(XP,YP,X,Y,XN,YN)  :- 
up_tum(UL),YP>Y,Y=YN, 
member((XU,YU),UL), 
YU>Y,X>=XU, 

add_statement(bump(XU,YU,X,YU)), 
delete((XU,YU),UL,ULl), 
add_statement(up_tum(UL  1 )), 
del_statement(up_turn(UL)), 
write_to_screen(XU,YU,X,YU),stars, 
repeat,check_duplicate(XP,YP,X,Y,XN,YN) 
write_to_file(XU,YU,X,YU),!. 


/*  rtn_to_horiz  is  below  latest*/ 
/*    up_tum,  so  create  the     */ 
/*   rtn_to_horizontal  coord.  */ 


create_upturn(XP,YP,X,Y,XN,YN)  :- 

Y  +  3<YP 

Y  =  YN, 

X  +  2  <=  XN, 

up_turn(UL), 

find_near_up_turn(UL,Y,FX,FY), 

COORD  =..  [cp,FX,FY], 

points(L), 

locate_coords(COORD,L,LL), 

find_crossing(LL,Y,XCOORD), 

write_to_file(XCOORD,Y,X,Y), 

write_to_screen(XCOORD,Y,X,Y). 


/*  rtn_to_horiz  is  below  latest  */ 
/*    up_turn,  so  create         */ 
/*    the  upturn.  */ 


check_duplicate(XP,YP,X,Y,XN,YN):- 

below_upturn(XP,YP,X,Y,XN,YN),!. 
check_duplicate(XP,YP,X,Y,XN,YN). 

find_near_up_turn([(FX,FY)IRL],Y,FX,FY)  :- 

FY<Y,  !. 
find_near_up_turn([AIRL],Y,FX,FY)  :- 

find_near_up_turn(RL,Y,FX,FY). 


/*  Make  sure  the  rtn_to_horiz  */ 
/*   doesn't  include  multiple  */ 
/*    bumps.  */ 


38 


fmd_crossing([A,BILL],Y,XCOORD):-  /*  Returns  vertical  coord     */ 

A  =..  [cp,AX,AY],  AY  <  Y,  /*    where  up.turns  Y  coords  */ 

B  =..  [cp,XCOORD,BY],  BY  >=  Y,!.  /*    was  crossed.  */ 

find_crossing([A,BILL],Y,XCOORD)  :- 

find_crossing([BILL],Y,XCOORD). 

locate_coords(COORD,[AIL],[AIL]) :-  /*  Find  a  coordinate  from  list*/ 

first_one(COORD,[AIL]),!. 
locate_coords(COORD,[AIL],LL)  :- 

locate_coords(COORD,L,LL). 

first_one(X,[XIL]). 

write_to_screen(X,Y,Xl,Yl)  :- 
stars, 

writeCbump  at  '),write(X),write('  '),write(Y),write('  '), 
write(Xl),write('  '),write(Yl),nl. 

write_to_file(X,Y,Xl,Yl)  :- 
tell(l, "bumps. pro"), 
write('b('), 
write(X),write(7), 
write(Y),write(7), 
write(Xl),  write(','), 
write(Yl),write('),'),nl, 
told(l). 

close_bump_file  :- 

up_turn(UL), 

close_up_turns(UL), 

tell(l, "bumps. pro"), 

write('b(999,999,999,999)]).'),nl, 

writeCendmod  /*bumps.pro  */'), 

told(l). 

close_up_tums([]).  /*  Closes  any  up_turns  left  */ 

close_up_tums([(X,Y)IL]) :-  /*    in  the  list.  */ 

end_pt([XE,YE]),write_to_file(X,Y,XE,Y), 

close_up_tums(L). 

endmod  /*  process  */  . 
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UTILITIES.PRO 

This  program  contains  utility  predicates  and  must  be  loaded  as  a  module  with 
PROCESS.PRO  and  BUMPID.PRO. 

module  utilities. 

export  (append  /  3,  add_item  /  3,  delete  /  3,  reverse  /  2,  singlemember  /  2, 

member  /  2,  stars  /  0,  first  /  2,  index  1  /  5, 
open  /  2,  tell  /  2,  told  /  1 ,  get_fm_file  /  2,  read_file  /  0). 

/*$eject*/ 

body. 

dynamic  (counter/  1). 

/*  These  file  were  written  by  Dave  Hutson,  originally  in  C-prolog.        */ 

/*  It  has  been  converted  to  run  on  M-prolog  by  Robert  Powell.  Other  utility  predicates  */ 

/*  have  also  been  added  as  needed.  */ 

abs(X,X)  :-  X  >=  0. 
abs(X,Y)  :-  Y  is  0  -  X. 

first([XIL],X). 

last([X],X). 

last([XIL],Y) :-  last(L,Y). 

member(X,[XIL]) :-  !. 
member(X,[YIL]) :-  member(X,L). 

delete(X,[],[]). 

delete(X,[XIL],M)  :-  !,  delete(X,L,M). 
delete(X,[YIL],[YIM]) :-  delete(X,L,M). 

append([],L,L). 

append([XIL],L2,[XIL3])  :-  append(L,L2,L3). 

add_item(X,L,[XIL]):-  !. 

reverse(L,R)  :-  reverse2(L,[],R). 

reverse2([],L,L) :-  !. 

reversc2([XIL],R,S)  :-  reverse2(L,[XIR],S). 

/*  writes  60  stars  */ 

stars  :- 

writeC***********  ******************************************  ******") 
nl,!. 
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I*         These  predicates  are  from  Professor  Rowe's  book.  */ 

I*  Iterates  repeatedly  forward  through  a  predicate  until  condition     */ 
/*  "done",  defined  by  the  user,  is  satisfied.  */ 

iterate(PRED)  :- 

repeat,  iterate2(PRED),  done  . 

iterate2(PRED)  :- 

evaluate(PRED), !  . 
iterate2(PRED) . 

I*         Written  by  Mike  Bizer  for  PROCESS.PRO  */ 

/*  indexes  into  a  list,  returning  the  previous(PP),item(P),  and  next(PN)  */ 

indexl([PP,P,PNIL],N,PP,P,PN)  :- 

N  =  0,!. 
indexl([AIL],N,PP,P,PN)  :- 
NCisN-1, 
index  1(L,NC,PP,P,PN),!. 

index(X,[XIL],l)  :- 
I 

index(X,[YIL],N)  :- 

index(X,L,NMl),  N  is  NM1  +  1  . 

/*  These  utilities  were  wrtten  by  Roben  Powell.  */ 

open(CH,OUTFILE):- 

set_channel(outfile(CH),[name=OUTFTLE,mode=create]), 

set_output(outfile(CH)), 

told(CH). 

tell(CH,OUTFILE):- 

set_channel(outfile(CH),[name=OUTFILE,mode=append]), 
set_output(outfile(CH)). 

told(CH)  :- 

close_output(outfile(CH)). 

get_fm_file(CH,INFILENAME)  :- 

set_channel(infile(CH),name=INFILENAME), 

set_input(infile(CH)), 

read_file, 

fail. 
get_fm_file(CH,INnLENAME)  :- 

close_input(infile(CH)). 


read_file  :- 

read_token(file  end),!. 
read_file  :- 

read(X), 

add_statement(X),   /*  bottom  */ 

read_file. 

endmod  /*  utilities  */  . 
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BUMPID.PRO 

This  program  identifies  the  bumps  found  by  PROCESS. PRO.    It  uses  the 
COORDS.PRO  and  BUMPS.PRO  data  files. 

module  bumpid. 

/*$eject*/ 
body. 

dynamic  (index_counterl  / 1). 

dynamic  (the_bump  /  5).         /*  bump  fact  to  be  processed  */ 

dynamic  (low/  1). 

dynamic  (left/  1). 

dynamic  (high/  1). 

dynamic  (right/  1). 

dynamic  (bump_coords  /  1)     /*  coordinates  making  bump  */ 

dynamic  (bump_descr  /  3).      /*  bump  &  description  list  */ 

dynamic  (box  /  8).  /*  box  description         */ 

dynamic  (finish_bump  /  1). 

dynamic  (bumps_to_process  /  1 ). 

dynamic  (shipjength  /  1). 

go  :-    state(cpu_time,TI), 

add_statement(index_counterl  (0)), 
set_state(global_stack,30000), 
set_state(main_stack,  10000), 
system(garbage_collection), 
system(compress_stacks), 
set_up_bumpid,iterate(gol ), 
state(cpu_time,TF),  TOTAL_TIME  is  TF  -  TI, 
stars,write("cpu  time  =  "),write(TOTAL_TIME), 
nl. 

gol  :-  nl, 

pbump(B  1  ),compute_shape(B  1  ),check_number_id(B  1 ), 

tell(l, "objects"), 

bump_descr(B  1  ,L,0), 

write(Bl),nl, 

writeC    BUMP_DESCR -->'), 

write(L),nl, 

writeC    BUMP_ID --> '), 

write(0),nl, 

told(l), 

clean_corners,index_counter  1  (X), 

X2  is  X  +  1,  del_statement(index_counterl(X)), 

add_statement(index_counter  1  (X2)),nl, 

check_statement_table(B  1 ), 

stars. 
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check_statement_table(Bl) :-  state(statement_table,[U,S]),U>300000EO, 
add_statement(finish_bump(b(999,999,999,999))),stars, 
write(" Statement  table  full  after  bump  "),write(Bl),stars,!. 

check_statement_table(B  1). 

clean_corners:- 

low(  Y 1 7),del_statement(low(Y  17)), 

high(Y18),del_statement(high(Y18)), 

right(X  1 8),del_statement(right(X  1 8)), 

left(X17),del_statement(left(X17)), 

bump_coords(XYZ),length(XYZ,N),write("  #  segments  ==>  "), 

write(N),nl, 

del_statement(bump_coords(XYZ)), 

add_statement(bump_coords([])). 

done  :-  finish_bump(B),bconvert(B,X,Y,Xl,Yl), 
X=999,  XI  =999. 

display_stats  :-  nl,stars, 

state(main_stack,[U,C]), 

write("main_stack  used  =  "),write(U),nl, 

state(statement_table,[U  1  ,C  1  ]), 

write ("statement_table  used  ="),write(Ul),write("  "),write(Cl), 

nl. 

pbump(Bl) :-  get_two_bumps(Bl,B2),asserta(finish_bump(B2)), 
convertit(B  1  ,X,Y,X  1 ,  Y 1  ),convertit(B2,X2,  Y2,X3,  Y3), 
add_statement(the_bump(B  1  ,X,  Y,X  1 ,  Y 1 )), 
add_statement(bump_descr(Bl, [],[])), 
create_bump_coords , 


set_up_bumpid  :- 

open(l, "objects"), 

tell(l,"objects"),nl, 

write("OBJECTS  from  BUMPID"),nl, 

told(l), 

bumps(B),add_statement(bumps_to_process(B)), 

start_pt([X,Y]),end_pt([Xl,Yl]), 

Z  is  XI  -  X,  add_statement(ship_length(Z)),nl, 

add_statement(finish_bump([] )), 

add_statement(bump_coords  ([])),!. 

get_two_bumps(Bl,B2)  :- 

bumps_to_process([Bl,B2IB]),  write(Bl),  write(' '),  write(B2),  nl, 
del_statement(bumps_to_process([B  1  ,B2IB])), 
add_statement(bumps_to_process([B2IB])),!. 
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convertit(B,X,Y,Xl,Yl)  :- 

B=..[b,X,Y,Xl,Yl]  . 
cpconven(B,X,Y)  :- 

B=..[cp,X,Y]. 
bconvert(B,X,Y,Xl,Yl)  :- 

B=..[b,X,Y,Xl,Yl]. 

create_bump_coords  :-  the_bump(B,X,Y,Xl,Yl), 

points(L),find_start(X,Y,L,LR), 

create_list(X,  Y,X  1 ,  Y 1  ,LR), 

del_statement(the_bump(B  ,X,  Y,X  1 ,  Y 1 )), ! . 
create_bump_coords  :-  the_bump(B>X>Y,Xl,Yl), 

points(L),rcverse(L,RL), 

find_rev_start(X  1 ,  Y 1  ,RL,LR), 

reverse(LR,RRL), 

create_list(X,  Y,X  1 ,  Y 1  ,RRL), 

bump_coords([FIC]), 

cpconvert(F,XB,YB), 

cpconven(D,XB,Yl),append(C,[D],CC), 

add_statement(bump_coords(CC)), 

del_statement(bump_coords([FIC])), 

del_statement(the_bump(B,X,Y,Xl,Yl)), 

I 

find_start(X,Y,[AILR],LR)  :- 

A  =..[cp,Xp,Yp],X=Xp,Y=Yp, 

bump_coords(LL),append([A],LL,LLl), 

add_statement(bump_coords(LL  1  )),del_statement(bump_coords(LL)), ! . 
find_stan(X,Y,[AILR],LR)  :- 

cpconvert(A,Xp,Yp),X=Xp,Y<Yp,cpconvert(B,X,Y),not(member(B,[AILR])), 

bump_coords(LL), 

append([B],LL,LLl), 

add_statement(bump_coords(LLl)),del_siatement(bump_coords(LL)),!. 
find_stan(X,Y,[AILL],LR)  :-  find_start(X,Y,LL,LR),!. 

find_rev_start(X,Y,[AILR],LR)  :- 

A  =..[cp,Xp,Yp],X=Xp,Y=YP, 

bump_coords(LL),append(LL,[  A]  ,LL1 ), 

add_statement(bump_coords(LLl)),del_statement(bump_coords(LL)),!. 
find_rev_stan(X,Y,[AILL],LR)  :-  find_rev_start(X,Y,LL,LR),!. 

create_list(X,Y,Xl,Yl,[]):- 

write("  POSSIBLE  ERROR,  ran  out  of  coordinate  points  while"),nl, 

write("      processing  bump  "),write(X),write("  "),write(Y), 

write("    "),write(Xl),write("  "),write(Yl),nl,!. 
create.listCX.Y^Xl.Yl^AIL])  :- 

A=. .[cp,Xf,Yf],Yf  <  Yl,!. 
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create_list(X,Y,Xl,Yl,[AIL])  :- 

A=..[cp,Xf,Yf], 

XI  =  Xf,  Yl  >=  Yf,  bump_coords(LL),append(LL,[A],LU), 

add_statement(bump_coords(LLl)),del_statement(bump_C(X>rds(LL)),!. 
create_list(X,Y,Xl,Yl,[AIL])  :- 

/*  peeling  off  upper  bumps  */ 

bumps(B), 

cpconvert(A,XT,YT), 

bconvert(TEST,XT,YT,XX,YY), 

member(TEST,B ), 

locate_coords(TEST,B,[BUMP_COORDSILLL]), 

bconvert(BUMP_COORDS^A,YA,XF,YF), 

bconvert(WORKING_BUMP,X,  Y,X  1 ,  Y 1 ), 

not(member(BUMP_COORDS,[WORKING_BUMP])),  /*ensure  not  same 

bump*/ 

cpconvert(C,XF,YF), 

member(C,[AIL]), 

locate_coords(C,[AIL],[ZIRL]), 

write("  peeling  off  bumpl  M),write(BUMP_COORDS),nl, 

create_list(X ,  Y,X  1 ,  Y 1  ,RL), ! . 

create_list(X,Y,Xl,Yl,[AIL]) :-  /*  don't  compare  with  bump  using  added*/ 

/*  peeling  off  upper  bumps  */ 
bumps(B),  /*  with  created  endings    */ 

cpconvert(A,XT,YT), 
bconvert(TEST,XT,YT>XX,YY), 
member(TEST,B), 

locate_coords(TEST,B,[BUMP_COORDSILLL]), 
bconvert(BUMP_COORDS,XA,YA,XF,YF), 
bconvert(  WORKING_BUMP,X,  Y,X  1 ,  Y 1 ), 
not(member(BUMP_COORDS,[WORKING_BUMP])),  /*ensure  not  same 

bump*/ 
cpconvert(C,XF,YF), 

locate_fake_coords(X  1 ,  Y 1  ,C,[  AIL],[ZIRL]), 
write("  peeling  off  bump2  "), 
write(BUMP_COORDS),nl, 
create_list(X,  Y,X  1 ,  Y 1  ,RL),! . 

create_list(X,Y,Xl,Yl,[A!L])  :- 

bump_coords(LL),append(LL,[A],LLl), 

del_statement(bump_coords(LL)),add_statement(bump_coords(LLl)), 

create_list(X,Y,Xl,Yl,L),!. 
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compute_shape(Bl)  :-  bump_coords([LLOWIL]), 
set_defaults([LLOWIL]), 
check_defaults(Bl), 
find_corners(L), 

make_box(XMIN,YMIN,XMAX,YMAX), 
write("#####  BOX  >  "), 
writeCXMIlVO.writeCYMnsO.writeC  '), 
writeCXM^/vmteCYMAXXwriteO '), 
write(XMAX),write(YMAX),write(* '), 
write(XMAX),write(YMIN),nl, 

describe_shape(B  1  ,[LLOWIL]^IN,YMm,XMAX,YMAX), 
identify_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX), 
I. 

make_box(XMIN,YMIN,XMAX,YMAX)  :- 
low(YMIN), 
left(XMIN), 
high(YMAX), 
right(XMAX), 
add_statement(box(XMIN,YMIN^MIN,YMAX,XMAX,YMAX,XMAX,YMIN) 

find_min(A,B,A)  :-  A  <=  B,!. 
find_min(A,B,B). 
find_max(A,B,A)  :-  A  >=  B,  !. 
find_max(A,B,B). 

set_defaults([LLOWIL]) :-  LLOW  =..  [cp,X,Y], 
add_statementGow(Y)), 
add_statement(left(X)), 
get_last([LLOWIL]),!. 

get_last([AI[]])  :-  A  «..[cp,X,Y], 

add_statement(right(X)), 
add_statement(high(Y)),!. 
get_last([AIL]) :-  get_last(L). 

find_corners([])  :-  !. 
find_corners([AIL])  :- 

A  =. .  [cp,XA,YA],high(Y), 

Y<YA, 

add_statement(high(YA)), 

del_statement(high(Y)),find_comers(L),!. 


find_comers([AIL])  :- 

A  =  .  [cp,XA,YA],left(X), 
XA<X, 

add_statementOeft(XA)), 
del_statement(left(X)),find_corners(L), ! . 
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find_corners([AIL])  :- 

A  =..  [cp,XA,YA],right(X), 

XA>X, 

add_statement(right(XA)), 

del_statement(right(X)),find_corners(L),! 
find_corners([  AIL]  ):- 

find_corners(L),!. 

check_defaults(B  1 )  :- 

bcon  vert(B  1  ,X,  Y.XX,  Y  Y), 

check_xmin(X), 

check_ymin(Y), 

check_xmax(XX), 

check_ymax(YY),!. 

check_xmin(X)  :- 

left(XT), 

XT  <=  X,  !. 
check_xmin(X)  :- 

left(XT), 

change_left_default(X,XT),! . 
check_ymin(Y)  :- 

low(YT), 

YT<=Y,  !. 
check_ymin(Y)  :- 

low(YT), 

change_low_default(Y,  YT), ! . 
check_xmax(XX)  :- 

right(XT), 

XT>=XX,  !. 
check_xmax(XX)  :- 

right(XT), 

change_right_default(XX,XT),!. 
check_ymax(YY)  :- 

high(YT), 

YT>=YY,  !. 
check_ymax(YY)  :- 

high(YT), 

change_high_default(YY,  YT),! . 

change_left_default(NEW,OLD)  :- 

add_statementGeft(NEW)), 

del_statement0eft(OLD)), ! . 
change_low_default(NEW,OLD)  :- 

add_statementGo  w(NEW)) , 

del_statementGow(OLD)), ! . 
change_right_default(NEW,OLD)  :- 

add_statement(right(NEW)), 

del_statement(right(OLD)),! . 
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change_high_default(NEW,OLD)  :- 
add_statement(high(NEW)), 
del_statement(high(OLD)),!. 

/*****************  LOCATION  descriptors  to  bump.descr  ***************/ 

describc_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX)  :- 
stan_pt([SX,SY]), 
ship_length(X), 

Y  is  (XMAX  +  XMIN)  /  2, 

Y  -  SX  <  15, 
add_descr(B  1  ,bow), 
describe_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX) ,  !. 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX):- 
start_pt([SX,SY]), 
ship_length(X), 

Y  is  (XMAX  +  XMIN)  /  2, 

Y  <  SX  +  (X  /  3), 

Y  -  SX  >=  15, 
add_descr(B  1  .forward), 
describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) ,  !. 

describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
stan_pt([SX,SY]), 
end_pt([EX,EY]), 
ship_length(X), 

Y  is  (XMAX  +  XMIN)  /  2, 

Y  >  SX  +  (2  *  X  /  3), 
EX  -  Y  >=  15, 
add_descr(Bl,aft), 
describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) ,  !. 

describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
start_pt([SX,SY]), 
ship_length(X), 

Y  is  (XMAX  +  XMIN)  /  2, 

Y  <  SX  +  (2  *  X  /  3), 

Y  >  SX  +  (X  /  3), 
add_descr(B  1  ,mid_ship), 
dcscribe_shape(Bl,[LLOWIL],XMlN,YMIN,XMAX,YMAX) ,  !. 

describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
end_pt([SX,SY]), 
ship_length(X), 

Y  is  (XMAX  +  XMIN)/ 2, 
SX  -  Y  <  15, 
add_descr(B  1  ,stern), 
describe_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX) ,  !. 
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/**%***************  acj^  sjze  descriptors  *********************/ 

describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

ship_length(X), 

OS  is  XMAX-XMIN, 

OS/X>=0.10E0, 

add_descr(B  1  ,extra_large), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX) , !. 
describe_shape(Bl,[LLOWIL],XMIN,YMIN)XMAX,YMAX)  :- 

ship_length(X), 

OS  is  XMAX-XMIN, 

OS/X  >=  0.03E0, 

OS/X<0.10E0, 

add_descr(B  1  ,large), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) , !. 
describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

ship_length(X), 

OS  is  XMAX-XMIN, 

OS/X  <  0.03E0, 

OS/X  >  0.0 15E0, 

add_descr(B  1  ,medium), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) ,  !. 
describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

ship_length(X), 

OS  is  XMAX-XMIN, 

OS/X  <  0.0 15E0, 

add_descr(B  1  ,small), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX) , !. 

/*****************  ^j  CURVINESS  to  bump_descr  ***************/ 

describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

bump_coords(XYZ),length(XYZ,N),N>50, 

add_descr(B  1  ,high_curviness), 

descrite_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) , !. 
describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

bump_coords(XYZ),length(XYZ,N),N<15, 

add_descr(B  1  ,no_curviness), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) ,  !. 

/*****************  add  descriptors  to  bump_descr  ***************/ 

describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

in_range(XMAX-XMIN,YMAX-YMIN), 

add_descr(B  1 ,  square), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) ,  !. 
describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

XMAX-XMIN  >  YMAX-YMIN  +  5, 

add_descr(B  1  ,flat_rectangle), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) , !. 
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describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

XMAX-XMIN  +  5  <  YMAX-YMIN, 

add_descr(B  l,tall_rectangle), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX) ,  !. 
describe_shape(Bl,[LLOWIL],XMIN,YMINfXMAX,YMAX):- 

XMAX-XMIN  <=  5,  (XMAX-XMIN)/(YMAX-YMIN)<=0.5E0, 

add_descr(B  l,pole), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX) , !. 
describe_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 

(XMAX-XMIN)/(YMAX-YMIN)<=0. 17E0, 

bump_coords(XYZ),length(XYZ,N),N<30, 

add_descr(Bl,pole), 

describe_shape(Bl,[LLOWIL],XMIN,YMIN>XMAX,YMAX) , !. 
describe_shape(Bl,[LLOWIL]^MIN,YMIN,XMAX,YMAX) :-  !. 

add_descr(Bl,X)  :- 

bump_descr(B  1  ,DL,0),not(member(X,DL)), 
append([X],DL,LL),add_statement(bump_descr(Bl,LL,0)), 
del_statement(bump_descr(B  1  ,DL.O)). 


/***********************  IDENTIFY  SHAPE  *****; 

identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descr(B  1  ,DL,6),member(pole,DL), 
not(member(bow,DL)),not(member(stern,DL)), 
add_id(B  1 , antenna), 
identify_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX),!. 

identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descT(B  1  ,DL,0),member(pole,DL), 
not(member(mid_ship,DL)),not(member(forward,DL)), 
not(member(aft,DL)), 
add_id(B  1  ,mast/support), 
identify_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX),!. 

identify_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX):- 
bump_descr(B  1  ,DL,0),member(flat_rectangle,DL), 
not(member(small,DL)),not(member(extra_large,DL)), 
meraber(high_curviness,DL), 
not(member(mid_ship,DL)), 
add_id(B  1  ,gun-turret/weapon-system), 
identify_shape(Bl,[LLOWIL],XMIN,YMIN^MAX,YMAX),!. 


%***********************/ 
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identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descr(B  1  ,DL,0),member(flat_rectangle,DL), 
not(member(mid_ship,DL)),  not(member(small,DL)), 
not(member(extra_large,DL)), 
not(member(high_curviness,DL)), 
not(member(no_curviness,DL)), 
add_id(B  1  ,  weapon-system), 
identify_shape(Bl,[LLOWIL],XMIN,YMINrXMAX,YMAX),!. 

identify_shape(B  1  ,[LLOWIL],XMIN,YMm,XMAX,YMAX)  :- 
bump_descr(B  1  ,DL,0),member(flat_rectarigle,DL), 
member(forward,DL),  not(member(small,DL)), 
(member(no_curviness,DL)), 
add_id(B  1  superstructure), 
identify_shape(Bl,[LLOWIL]^MIN,YMIN^MAX,YMAX),!. 

identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descr(Bl,DL,0), 
member(mid_ship,DL), 

not(member(pole,DL)),not(member(small,DL)), 
not(member(tall_rectangle,DL)), 
(XMAX-XMIN)/(YMAX-YMIN)  <=  2.0E0, 
(XMAX-XMIN)/(YMAX-YMIN)  >=  0.5E0, 
add_id(Bl  .radar), 
identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX),! . 

identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descr(B  1  ,DL,0), 

member(mid_ship,DL),  not(member(no_curviness,DL)), 
not(member(pole,DL)),not(member(square,DL)), 
not(member(tall_rectangle,DL)), 
(XMAX-XMIN)/(YMAX-YMIN)  <=  2.0E0, 
(XMAX-XMIN)/(YMAX-YMIN)  >=  0.5E0, 
add_id(Bl,radar), 
identify_shape(Bl,[LLOWIL],XMIN,YMIN,XMAX,YMAX),!. 

identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descr(B  1  ,DL,0),member(tall_rectangle,DL), 
not(member(pole,DL)),not(memberGarge,DL)), 
(XMAX-XMIN)/(YMAX-YMIN)  <=  2.0E0, 
add_id(B  1  ,mast/support), 
identify_shape(B  1  ,[LLOWIL]^CMIN,YMIN^MAX,YMAX),! . 

identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
bump_descr(B  1  ,DL,0),not(member(pole,DL)), 
not(member(small,DL)),not(member(high_curviness,DL)), 
not(member(forward,DL)), 
add_id(B  1  superstructure), 
identify_shape(B  1  .[LLOWIL]  ,XMIN,  YMIN,XMAX,  YMAX), ! . 
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identify_shape(B  1  ,[LLOWIL],XMIN,YMIN,XMAX,YMAX)  :- 
check_if_empty  (B 1 ),! . 

add_id(Bl,X):- 

bump_descr(B  1  ,DL,0),not(member(X,0)), 

append([X]  AOL),add_statement(bump_descr(B  1  ,DL,OL)), 

del_statement(bump_descr(B  1  ,DL,0)). 

/************  change  list  to  OR  if  multiple  identifications  made  ****/ 

check_number_id(Bl)  :- 

bump_descr(Bl,DL,0),length(0,N),write(N),N>l, 

write(0),nl,change_to_or(0,or,00), 

add_statement(bump_descr(B  1  ,DL,00)), 

del_statement(bump_descr(B  1  ,DL,0)),!. 
check_number_id(B  1 ) :- !. 

change_to_or([A,BIL],OR,LL)  :- 

append([OR],[A],Y),append([B],Y,LL),!. 

check_if_empty(Bl)  :- 

bump_descr(Bl,DL,[]), 

add_id(B  1  .unknown), ! . 
check_if_empty  (B 1 ). 

locate_coords(COORD,[AIL],[AIL])  :- 

first_one(COORD,[AIL]),!. 
locate_coords(COORD,[AIL],[AIL])  :- 

cpconvert(COORD,X,Y), 

cpconvert(A,XT,YT), 

X=XT,  Y  >  YT,!. 
locate_coords(COORD,[AIL],LL)  :- 

locate_coords(COORD,L,LL). 

locate_fake_coords(Xl,Yl,COORD,[AIL],[AIL]):- 

cpconven(COORD,XF,YF),cpconvert(A,XT,YT), 

XT=XF,YT<YF, 

bump_coords(BL),  append(BL,[COORD],BLL),append(BLL,[A],LL), 

add_statement(bump_coords(LL)), 

del_statement(bump_coords(BL)), 

;# 

locate_fake_coords(X  1  ,Y1  ,COORD,[AIL],LL)  :- 

locate_fake_coords(X  1 ,  Y 1  ,COORD,L,LL). 

first_one(X,[XIL]). 
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/*************************  utilities  ********************************/ 
in_range(VARl,VAR2):- 

VARl<=VAR2  +  3, 

VAR1  >=  VAR2  -  3. 

endmod  /*  bumpid  */  . 
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APPENDIX  C  ■  OUTPUT  FROM  GARCIA  (LAMPS) 


GAR-OBJ  OUTPUT 

This  is  an  example  of  output  from  the  GARCIA  (LAMPS)  silhouette,  showing  the 
bump  locations,  description,  and  identification. 

OBJECTS  from  BUMPID 
b(341,484,361,484) 

BUMP_DESCR-->[flat_rectangle,high_curviness,large,forward] 

BUMP_ID  -->  [gun-turret/weapon- system] 
b(340,476,364,476) 

BUMP_DESCR  — >[flat_rectangle,no_curviness,large,forward] 

BUMP_ID  -->  [superstructure] 
b(407,485,415,485) 

BUMP_DESCR-->[flat_rectangle,large,forward] 

BUMP_ID  -->  [weapon- system] 
0(403,478,416,478) 

BUMP_DESCR  -->[flat_rectangle,no_curviness,medium,forward] 

BUMP_ID  -->  [superstructure] 
0(457,537,459,537) 

BUMP_DESCR  -->[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMPJD  -->  [antenna] 
b(456,522,457,522) 

BUMP_DESCR  -->[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMP_ID  -->  [antenna] 
b(465,564,465,564) 

BUMP_DESCR  -->[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMP_ID  -->  [antenna] 
b(465,522,467,522) 

BUMP_DESCR-->[pole,tall_rectangle,small,mid_ship] 

BUMP_ID  -->  [antenna] 
0(489,526,499,526) 

BUMP_DESCR  -->[square,no_curviness,small,mid_ship] 

BUMPJD  -->  [unknown] 
0(533,543,536,543) 

BUMP_DESCR  -->[flat_rectangle,high_curviness,large,mid_ship] 

BUMP.ID  -->  [radar] 
0(536,571,541,571) 

BUMP_DESCR  ">[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMP_ID  -->  [antenna] 
b(559,618,561,618) 

BUMP_DESCR  -->[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMP_ID  -->  [antenna] 
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b(565,618,569,618) 

BUMP_DESCR~>[square,no_curviness,small,mid_ship] 

BUMP_ID  -->  [unknown] 
b(560,60 1,56 1,601) 

BUMP_DESCR-->[tall_rectangle,small,mid_ship] 

BUMP_ID  -->  [mast/support] 
b(549,592,563,592) 

BUMP_DES CR  -->[flat_rectangle,large,mid_ship] 

BUMP_ID  ->  [superstructure] 
b(529,539,553,539) 

BUMP_DES CR -->[flat_rectangle,no_curviness,large,mid_ship] 

BUMP_ID  -->  [superstructure] 
b(552,544,553,544) 

BUMP_DESCR-->[tall_rectangle,medium,mid_ship] 

BUMP_ID  -->  [mast/support,or,superstructure] 
b(551,571,553,571) 

BUMP_DESCR  -->[square,medium,mid_ship] 

BUMP_ID  -->  [radar,or,superstructure] 
b(527,528,554,528) 

BUMP_DESCR-->[flat_rectangle,large,mid_ship] 

BUMP_ID  -->  [radar,or,superstructure] 
0(485,514,559,514) 

BUMP_DESCR-->[flat_rectangle,extra_large,mid_ship] 

BUMP_ID  -->  [superstrucnire] 
0(587,507,589,507) 

BUMP_DESCR-->[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMP_ID  -->  [antenna] 
b(602,502,617,502) 

BUMP_DESCR-->[flat_rectangle,no_curviness,medium,mid_ship] 

BUMP_ID  -->  [superstrucnire] 
0(631,516,633,516) 

BUMP_DESCR-->[pole,tall_rectangle,small,mid_ship] 

BUMP_ID  -->  [antenna] 
b(626,502,629,502) 

BUMP_DESCR-->[tall_rectangle,small,mid_ship] 

BUMP_ID  -->  [mast/support] 
b(65 1,494,653,494) 

BUMP_DESCR  ->[pole,tall_rectangle,no_curviness,small,mid_ship] 

BUMP_ID  -->  [antenna] 
b(675,496,699,496) 

BUMP_DESCR-->[flat_rectangle,high_curviness,large,aft] 

BUMP_ID  — >  [gun-turret/weapon-system] 
0(737,515,739,515) 

BUMP_DESCR  -->[pole,tall_rectangle,no_curviness,small,aft] 

BUMPJD  -->  [antenna] 
b(745,520,747,520) 

BUMP_DESCR-->[pole,tall_rectangle,no_curviness,small,aft] 

BUMP_ID  -->  [antenna] 
b(745,515,755,515) 

BUMP_DESCR  -->[tall_rectangle,medium,aft] 

BUMP_ID  -->  [mast/support,or,superstructure] 
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b(788,517,798,517) 

BUMP_DESCR  -->[no_cuiviness,small,aft] 

BUMP_ID  -->  [unknown] 
b(763,515,802,515) 

BUMP_DESCR-->[flat_rectangle,no_curviness,large,aft] 

BUMP_ID  -->  [superstructure] 
b(745,514,803,514) 

BUMP.DESCR  -->[square,large,aft] 

BUMP_ID  -->  [superstructure] 
0(675,494,803,494) 

BUMP_DESCR-->[flat_rectangle,high_curviness,extra_large,aft] 

BUMP_ID  -->  [unknown] 
0(725,499,803,499) 

BUMP_DESCR-->[flat_rectangle,extra_large,aft] 

BUMP_ID  -->  [superstructure] 
0(932,504,933,504) 

BUMP_DESCR-->[pole,no_curviness,small,stern] 

BUMP_ID  -->  [mast/support] 
b(93 1,497,93 1,497) 

BUMP_DESCR  ->[pole,no_curviness,small,stem] 

BUMP_ID  -->  [mast/support] 
b(929,49 1,934,491) 

BUMP_DESCR  -->[pole,tall_rectangle,no_curviness,small,stem] 

BUMP_ID  -->  [mast/support] 
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